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Abstract: Nowadays, climate models rely on couplers. Each complete cli- 
mate model is broken into different sub-models (oceanic, atmospheric,...), 
each one working on a different grid. The coupler brings these models to- 
gether and interpolates the physical quantities between the grids. However, 
neither the coupler nor sometimes the sub-models themselves know precisely 
the grid cell edges. They only know the grid cell corners (vertices) and the 
true grid cell areas. Thus, the coupler has to make assumptions about the 
grid cell edges in order to compute the grid cell intersections. For first-order 
schemes, the most straightforward way to interpolate scalar quantities is to 
directly use these approximate grid cell intersections, that don' take the true 
grid cell areas into account. It is the method used in the "conservative" re- 
gridding option implemented in the widely used spherical interpolation pack- 
age SCRIP. We show that is doesn't preserve integrals in the general case, 
whether the coupler using the SCRIP-generated weights transmits directly 
the intensive quantity F^j, or its extensive counterpart, FgiA^- (that is the 
same quantity multiplied by the true area of the individual source grid cell). 
We show how to modify the interpolation scheme to preserve integrals in the 
general case. 

Key-words: conservative regridding, coupler, SCRIP, Lagrange multiplier, 
climate 
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1 Introduction 

SCRIP is a library widely used by different couplers between oceanic and 
atmospheric models [3^, '5]. In particular, it implements an option, the "con- 
servative" option, that should let the coupler preserve integrals of a quantity 
to be mapped between the two grids. 

SCRIP permits to map and interpolate a quantity F between two numer- 
ical grids in spherical coordinates (Fig. [T]). The only information given to 
SCRIP are the vertices of each source and destination grid cell. SCRIP, not 
knowing the true area of each grid cell, has to estimate these areas based on 
the coordinates of the vertices. As such estimated areas can differ from the 
true ones, the conservative option of SCRIP is mathematically wrong in the 
general case. 

Recently, different reports in the climate community have addressed this 
issue in U]. K.E.Taylor also suggested that the SCRIP library could be 
conservative in the general case, provided that the coupler multiplies the 
quantity to remap by the true area of the grid cell. |lj 

We will show under which conditions SCRIP is indeed conservative. We 
will give the minimum modifications that should be brought to make it con- 
servative in the general case. 

Apart from the use of the word 'SCRIP', which is the name of the in- 
terpolation library used in climate, all the mathematics of this report are 
general. It applies whenever we interpolate fields between two grids that 
have unknown grid cell edges. 

2 Scope of the report 

We want to know how the integral of a quantity can be preserved when 
interpolating between two grids. By hypothesis, the source map "contains" 
the destination map (see So, we don't analyze the issues linked to non- 
matching sea-land mask. 

3 Notations 

The quantity to be regridded is called F. F is an intensive quantity, such as 
the temperature. F is mapped between a source grid S and a destination grid 
D (Fig. [T|. All the quantities linked to the source grid take the underscript 
s (^s), whereas all the quantities linked to the destination grid take the 
underscript d (X^). The underscript i (Xj) always means "a given cell of the 
source grid", whereas the underscript j {Xj) always means "a given cell of 
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Figure 1: Scheme of the regridding : we interpolate from a source grid (un- 
derscript s) to a destination grid (underscript d). We note Asi, the area of 
the i*'* cell of the source grid. We note Adj, the area of the j*'^ cell of the 
destination grid. Prom Agi, A^j and the grid cell corners (vertices), we de- 
termine the interpolation weights fij. We treat first-order schemes, so the 
intensive quantity Fgi is constant over the i*'^ source grid cell. The inter- 
polated quantity F^j is constant over the j*'' destination grid cell. It reads 

^dj = fijFsi- 

the destination grid". The superscript T {X'^) stands for true, whereas the 
superscript F {X^) stands for false. The letter A stands for area. Ai fl A2 
stands for the area of the intersection of the two corresponding surfaces Si 
and 5*2. 
Thus, 

• ^4^ is the true area of the i*'^ cell in the source grid, that is the area 
used by the model on the source grid, 

• A^p the true area of the j*^ cell in the destination grid, 

• is the estimated area of the i*'^ cell in the source grid. The estimation 
is done and used inside SCRIP. We say that A^^ is the false area of the 

• A^-, the estimated (false) area, of the j*'* cell in the destination grid, 

• Fsi, the constant value of the field P in the i*'^ cell of the source grid, 

• Fdji the regridded value of F in the j*'^ cell of the destination grid. The 
coupler estimates F^j from F^,, 

• (j)si = X Fsi J integral of Fgi over the true area of the i*'* cell of the 
source map. It is an extensive variable. 
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• (f)dj = X Fdj, integral of F^j over the true area of the j* cell of the 
destination map. It is an extensive variable, 

From now on, we use 'SCRIP' as a shortcut of 'the conservative regridding 
option implemented in the SCRIP'. 



4 Definitions 

Definition 1 A regridding scheme calculates the sparse matrix that trans- 
forms Fsi into Fdj . The coefficients of the matrix are called fij : 



Definition 2 Alternatively, a regridding scheme can calculate the sparse ma- 
trix that transforms (f)si into ■ The coefficients of the matrix are called Pij : 



(t>dj = ^Pi34>si, (2) 



Pij derives from fij : 



Pij — fij X (3) 



SI 



Definition 3 The regridding of F is globally conservative if and only if 
J2i'Psi = J2j4^dj- From now on, we use 'conservative' in place of 'globally 
conservative '. 

Definition 4 The regridding of F from a source to a destination grid is 
locally conservative if and only if the flux on the destination grid cell is exactly 
equal to the sum of all the fluxes from the underlying source grid cells. If the 
intensive quantity Fgi is constant inside each source grid cell, this writes 
exactly 

<t>d, = Y}^siX{AlnAlj)]. (4) 

i 

This gives the two equivalent expressions for a locally conservative regridding 
scheme 

'Pdj = ^[4>si X T '^' \ (5) 

si 
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and 



\Ps^ X ^^], (6) 



which, given Eq. [7], writes 



jij at ■ y' ' 

eq:regridding 

Note that a locally conservative regridding scheme is also conservative 
(see the demonstration in Eq. 16). 

Definition 5 The regridding doesn't introduce fluctuations if and only if a 
constant field is transformed into the same constant field. This can also he 
expressed as : Fgi = F,'^i implies Fdj = F,'^i, or, by using Eq. [I] ; 



5 Implicit hypotheses 

In this report, we implicitly make the following assumptions over the source 
and the destination map : 

Hypothesis 1 Each true destination grid cell is included in a subset of true 
source grid cells; the true grid cells are pairwise disjoint : 

niAl = and = 

This implies that 

E^^^^^5 = 4 (9) 

i 

Hypothesis 2 Each false source grid cell is included in a subset of false 
destination grid cells; the false grid cells are pairwise disjoint : 

Hi/l^ = and = 

This implies that: 

E^SnA5 = /if, (10) 
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6 Analysis of current SCRIP 

The current version of SCRIP calculates /jj (Def. [T]) as : 

^SCRIP _ si I ' ^dj /. -, N 

Jij aF \^^) 

Note that, by definition, as for a locally conservative scheme (Def. |4]), 

i 

whereas, in general, 

Y^ffr'^^- (13) 

3 

Note also that, from Eq. 12 Def. |5] and under Hyp. |2| SCRIP regrid- 



ding never introduces fiuctuations. Under SCRIP, a constant field is always 
regridded into the same constant field. 

6.1 False and true grid areas are identical 

In this paragraph, we make the following hypothesis : 
Hypothesis 3 The false and the true grid cell areas are identical, that is 

= Al (14) 

and 

Adj = A'^j. (15) 
This means that SCRIP is locally conservative (see Eq. 11 and Def. |^. 
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j j 



= j:(J2f'r'Fs.)xAl fromEq.[T] 



(16) 



from Eq. [TT] and Eq. 14 15 



J2Fs^Al from Hyp. [T] 



i 



So, given Def. [3} if the areas are estimated correctly, SCRIP is conserva- 



tive. 



6.2 False and true grid areas are different 
6.2.1 The coupler using SCRIP transmits Fgi 

We now test if SCRIP is conservative when the false and the true grid areas 
are different : 



j i 3 

= E[^-x5^(A^nAj.)x^] from Eq. [11] (17) 



So, SCRIP is not conservative. 

6.2.2 The coupler transmits (psi = FsiA^^ instead of Fgi 

We now consider the case in which the coupler transmits the extensive quan- 
tities ipsi = FgiA'^- instead of Fgi. This means that we apply the following 
recipe : 

1. we don't modify the SCRIP code, 
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2. the source map sends (psi instead of Fsi to the destination map. 

3. in the destination map, we recover Fdj = 4f- by dividing the output 

of the matrix multiphcation with the true area of the destination grid 
cell 



Eq. [T] becomes 



which also writes 



^^.4 = E/™^-^- (19) 

i 

= j:if"r' X (20) 

We want to determine if such a scheme makes SCRIP conservative : 



(21) 



So, 



4» 



4» 



Z1'^^^ = 5Z'^- (22) 

j i 

Y.f"r' = ^ V. (23) 

j 

E^'^^l « (24) 



But, ^ . ".f ''^ 7^ 1 in the general case. So, transmitting (f)si in place of 



^3 A 

Fsi doesn't render SCRIP conservative. 



7 Modifications needed to render SCRIP con- 
servative 

We propose the modifications that need to be brought to SCRIP or to the 
use of SCRIP to make the regridding conservative in the general case. 
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7.1 The coupler transmits the intensive quantities Fsi 

We are looking for the new weights y^^CRiP.mt ^^isX would make SCRIP con- 
servative, when the coupler transmits the intensive quantities Fg^. 



From Eq. 17, SCRIP is conservative if and only if 



E4' 



SCRIP,int 



Ai 



(25) 



We know that 



1 (see Hyp. 



2). So, 



SCRIP,int 



is a solution of Eq. j25}^ 

We can write • ''° differently : 



(26) 



f 



SCRIP,int 



SCRIP ^ 1^ ^ -^J^ 
ij A^ A'T 

^si ^dj 
SCRIP ^li/^fi 



(27) 



So, if the coupler using SCRIP is not modified (the F^i are transmitted), 
SCRIP needs to know the true areas of the source and of the destination 
grid cells. We need to add a parameter (SCRIP_CON_SERVATIVE_- 



INTENSIVE) that tells SCRIP to compute fl 



SCRIP,int 



(Eq. 2Qh in place of 



This "intensive" modified SCRIP would always be conservative for inten- 
sive quantities. 



7.2 The coupler transmits the extensive quantities (j)si = 



F A^- 



Eq. 24 would always be true if A^^ was the denominator in place of A^-. 
We define 

A^ n A^ 

rSCRIP,ext _ ^si I ' ^dj 

J V 



A^ 



A^ 

^SCRIP 



P (2^ 
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instead of Fsi, then 



If the coupler transmits the extensive quantities 
SCRIP doesn't need to know the true grid areas. However, we need to add 
a parameter (for example SCRIP_CONSERVATIVE EXTENSIVE) that 
tells SCRIP to compute /^^-R^^'^^* (Eq. 



28) in place of ff^^^^ (Eq. 11). 



This "extensive" modified SCRIP woula always be conservative for exten- 
sive quantities. 



8 Discussion of the equations 

We have demonstrated that SCRIP as such is not conservative in the general 
case. Transmitting extensive instead of intensive quantities between the two 
grids doesn't render SCRIP conservative. We have modified the interpolation 
scheme to guarantee conservation. 

Having said this, we now discuss four points , the first ones being specific 
to the climate community, the last one having a more general scope : 

8.1 Effect of a global normalization 

When the grid cells are closer to the pole, the difference between true and 
false grid cell areas, as calculated by SCRIP, is larger. So, applying a global 
normalization would not resolve the imbalance between the pole and the 
equator introduced by the inaccurate conservative scheme. 

8.2 Working on different sphere surfaces 

Sometimes, the models to be coupled work on different sphere surfaces. A 
global normalization is enough to correct for such a mismatch. 

8.3 Why an accurate conservative regridding scheme? 

A small inaccuracy in the conservative regridding scheme is not random but 
it is both geographically localized and repeated at each time step of the 
simulation. This implies that it may alter significantly the result of the 
simulation and that the analysis of this accuracy is important. 

8.4 Numerical diffusion 

For the first order conservative regridding scheme (the one treated here), 
the quantity is distributed from the source grid cell to all destination grid 
cells intersecting it (see for a comparison between first- and second-order 
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remapping schemes the paper of Jones [2j). The quantity is distributed not 
only to the portion of the destination grid cell intersecting the source grid 
cell, but also to the whole surface of this destination grid cell. It acts as an 
numerical diffusion of the quantity to be mapped. It is unavoidable in the 
case of a first-order remapping scheme. 

When we modify the regridding scheme, we want to know if we increase 
this numerical diffusion. 

In the current SCRIP, this diffusion is not determined by the intersections 
between the true source grid cells and the true destination grid cells (as it 
would be in the ideal case), but by the intersections between the false source 
grid cells and the false destination grid cells. If we change the equations of 
SCRIP as proposed in this report, the intersections are between the same 
grid cells. So, we will not significantly modify the extent of the numerical 
diffusion, which is good. 



8.5 Numerical fluctuations 

In SCRIP, 

i 

So from Def. [5} the scheme doesn't introduce numerical fluctuations, that 
is a constant field is always transformed into the same constant field. But, 
Ylj ffj^^^^ 7^ 1, so, from Eq. 17, SCRIP is not conservative. 

In the new scheme (SCRIF1C0NSERVATIVE_INTENSIVE), 



jSCRIP,int ^ ^ 

j 

So, from Def. |4[ it is locally and globally conservative. However, 

i 

So, from Def. |5| it introduces numerical fluctuations. 



8.6 Being conservative and avoiding numerical fluctua- 
tions 

Mathematically, in a first-order scheme, the scheme avoids numerical fluctu- 
ations if and only if the sparse matrix of the interpolation weights [fij] is a 
left stochastic matrix. The scheme is locally conservative if and only if [fij] 
is a right stochastic matrix. In both cases, [ftj] is sparse. 



MDLS 



Conservative regridding 



13 



The scheme has both properties if and only if [fij] is a sparse doubly 
stochastic matrix, that is a matrix for which each row and column sums to 
1. We can state it as a constrained optimization problem and express it as 
Lagrange multipliers : we want to find the sparse doubly stochastic matrix 
[fij] that is the closest to the matrix [ffj"^^^]- The Lagrange function writes : 



~ -^ijifij fij ) 

= if 



(29) 



Bij is constant in this optimization problem. Its elements are set to force 
the sparsity of [fij]- Ideally, we want the final matrix [fij] to have the same 



sparsity as the matrix [/■^■' 



SCRIPl 

ij 



correctly distributed, the matrix [Bij] reads : 



B, 




f SCRIP ^ 

fSCRIP 
ij 







otherwise, we set more generally : 



So, if the degrees of freedom of [/j^*^^^^] are 



(30) 



Bij 


= 1 


if 


Bij 


< 1 


if 


Bij 


= 


if 



i close to j 
i far from j 

Then, [fij] is the solution of the linear system : 



dL 
dfij 
dL 
d\i 
dL 
dm 



= 
= 
= 
= 



if Bij 







if Bij ^ 
Vi 



(31) 



(32) 



The first-order interpolation using these interpolation weights is conser- 
vative and it avoids at the same time numerical fluctuations. As drawback, 
it can introduce a larger numerical diffusion. 



9 Conclusion 

The equations and the discussion presented above have a general scope. They 
apply whenever we want to make a first-order interpolation of a field between 
two grids, whose edges are not exactly known. 
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One way (the one used by SCRIP) is to estimate internally the edges and 
calculate the interpolation weights by directly computing the intersections of 
the estimated (false) areas of the source and of the destination grid cells. This 
method doesn't introduce numerical fluctuations but is not conservative. 

We have proposed two simple modifications to this algorithm : 

1. if intensive quantities are used by the coupler, SCRIP needs to know 
the true grid cell areas of the source and of the destination map, and to 



modify the algorithm to calculate the weights y^^CRiP,mt ^-^^^ 27)^ This 



"intensive" modified SCRIP would always be conservative for intensive 
quantities. 

if extensive quantities are used by the coupler, SCRIP needs to modify 



the algorithm to calculate the weights j^^^^^'"^^^ {^(l- 28). This "ex- 



tensive" modified SCRIP would always be conservativefor extensive 
quantities. 

Given the discussion of ^ we think that the climate community should apply 
these changes. [?] 



References 

[1] e-mail from K.E.Taylor. I'm not positive what SCRIP does exactly, but 
I would note that there is a relationship between the fraction of source 
areas intersecting target areas and the fraction of target areas intersecting 
source areas: f*_ji = f_ij * A'_si/A'_dj where A' are the areas as 
estimated by SCRIP. So, with my current understanding: What I have 
learned is that if I want to regrid and preserve the integral quantity, I 
can do this with any of the regridders we've been discussing by: 1. pre- 
processing my field by multiplying source grid cell values by the source 
grid cell areas (or fractions of areas in the case of partially masked cells). 
2. passing this pre-processed field to the regridder along with the vertices 
of the each source and target grid cell 3. post-processing the returned 
field by dividing the target grid cell values by the target grid cell areas 
(or fractions of areas in the case of partially masked cells). This would 
work I think (but haven't actually checked) when choosing the option of 
not passing grid cell areas to SCRIP. 

[2] Philip W. Jones. First- and Second-Order conservative remapping 
schemes for grids in spherical coordinates. Monthly Weather Review, 
127(9) :2204-2210, September 1999. 



MDLS 



Conservative regridding 



15 



[3] Philip W. Jones. A User's Guide for SCRIP: A Spherical Coordinate 



Remapping and Interpolation Package, http : //climate . laiil . gov/ 



Sof tware/SCRIP/SCRIPusers .pdf 2001. [version 1.4]. 



[4] Karl E. Taylor. Preserving Global Integrals When Regridding. Technical 
report, July 2010. 

[5] Li Yan and Sophie Valcke. Work report - implement the true area nor- 
malization in 0ASIS3. Technical report, 2011. 



RT 11° 0001 




Editeur 

MDLS - Maison de la Simulation USR 3441 - Bat. 565 - CEA Saclay - 91 191 Gif-sur-Yvette (France) 

http://www.maisondelasimulation.fr 



